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A generalization of the kinetic equation [l[ is proposed for explaining observed shapes of wind 
■ wave spectra. The approach allows to fix a critical uncertainty in modeling wind wave spectra using 

a trivial condition of equilibrium of nonlinear transfer and wave dissipation due to breaking. The 
first results of the generalized kinetic equation simulation are presented in order to show relevance 
of the proposed theory to the observed Phillips' spectrum E ~ ui~ for wind waves. 
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I. INTRODUCTION 



Wind waves in seas and lakes represent an extremely complicated natural phenomenon governed by a great number 
of physical mechanisms. Complete physical theories for many of these mechanisms are not proposed yet that makes 
particular asymptotic cases extremely valuable. Phillips Q considered a limiting case of saturation of a random wave 
field due to wave breaking. The dimensional consideration gave the famous Phillips energy spectrum E(w) — ag 2 u~ 5 . 
Phillips himself emphasized the idea of 'hard, saturated limit of the spectrum' where 'any excursion above this limit' 
*35 ' ' S ' re ^ evec ^ immediately by breaking' [see p. 506 of [3J. 

The conceptual passage from the idea of saturation to the concept of an equilibrium of random wave field brought 
Phillips himself Q to the model of spectral balance where all the mechanisms, wave input by wind, dissipation due 
to breaking and nonlinear transfer due to resonant wave-wave interactions compete on equal terms. This model lead 
to other fundamental exponent —4 of wind-wave spectral equilibrium range. The essential flaw of the Phillips theory 
I . of spectral equilibrium, in its own words [sect. 5 of 3], is a number of empirical parameters responsible for wind- wave 

■ spectral balance. All these parameters should be matched each other. Phillips formulated the corresponding criteria of 
. consistency of his model but pointed out the difficulties of this matching and deficiency of our knowledge of wind- wave 

■ physics. 

In contrast to the Phillips model Q of wind-wave spectra w -4 [as well as the semi-empirical model by Toba 0, 
that gives the same spectral slope] the weakly turbulent Kolmogorov-Zakharov cascade solutions [!, H| appear to 
' be 'self-sufficient' for explaining exponents —4 of direct and —11/3 of inverse cascades as inherent features of deep 
water wave spectra governed by constant spectral flux in an equilibrium range of wave scales. Thorough analysis of 
wind-wave balance shows that this equilibrium range where nonlinear resonant interactions are dominating does exist 
for typical conditions of wind seas and for wave scales we are really interested in 0-Q • It must be admitted that the 
scheme of wind-wave balance by Phillips Q remains the most recognized one in the wind-wave community in spite 
of a number of unresolved questions on this scheme correctness. The most critical questions concern mechanisms of 
wave input and dissipation. Today, our knowledge of these mechanisms is mostly empirical and their description is 
mostly based on quasi-linear parameterizations. Thus, one has a conceptual gap between this quasi-linear vision of 
spectra close to oj~ 4 and the 'utterly nonlinear' authentic Q spectrum w -5 . The present paper is aimed at filling this 
gap following the theory of wave turbulence [l(| and its recent extension [l[ . 

One of the point of critics of the weakly turbulent approach is restrictions on wave scales (or time of evolution) 
when the corresponding asymptotic expansions are formally valid. Generally, the basic assumptions of the theory fail 
to be valid at long time and/or very short or very long wave scales. At the same time, for particular solutions these 
basic assumptions keep their formal validity for infinitely long time in the whole range of wave scales. Newell and 
Zakharov [l[ proposed to consider the corresponding wave spectra as generalized Phillips spectra in wave turbulence. 

In this paper we use the extension of the wave turbulence theory [l| for explaining the authentic Phillips' spectrum 
lo~ 5 as a result of competition of wave-wave resonant interactions and inherently nonlinear dissipative processes. The 
new dissipation function can be found in physically consistent way as one corresponding to the generalized Phillips' 
spectrum. This new dissipation function has remarkable features. First, the balance of wave-wave resonances and 
nonlinear dissipation fixes the dissipation rate for the classic Phillips spectrum cj~ 5 . This surprising result is in 
dramatic contrast with conventional Phillips model where the dissipation rate is assumed to be arbitrary. Secondly, in 
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the generalized Phillips model of this paper the new dissipation absorbs completely the spectral energy flux directed 
to infinitely high wave numbers. Thus, the extended wave turbulence approach appears to be able to describe the 
whole range of wave scales and explain different exponents of the observed wave spectra [U EH • 

In § 2 we give basic equations and necessary comments to the concept of the generalized Phillips' spectra. A simple 
isotropic model of the generalized Phillips spectrum of deep water waves is discussed in § 3. Stationary and self-similar 
solutions for this model are analyzed as a basis for further comparison with results of simulations of the extended 
kinetic equation for deep water waves. Numerical simulation of wind-wave spectra within the kinetic equation with 
the new dissipation function are presented in § 4. Discussion and conclusions are given in § 5. 



II. GENERALIZED PHILLIPS' SPECTRUM OF DEEP WATER WAVES 

In this section we reproduce results of Newell & Zakharov [l[ for particular problem of statistical description of 
wind- wave field. We consider the Hasselmann kinetic equation for spatial spectrum of wave action of wind-driven 
waves in the following form 

dN k 



Of 



Vk^k V r iVk — S n l + Sin + Sdiss (1) 



Subscripts k, r for V are used for gradients in wavevector k and coordinate r spaces correspondingly. For Ak(i) 
and Uk the subscript k means dependence on wavevector. The term S n i in ([1]) is responsible for four-wave resonant 
interactions. Terms Si n and Sdiss describe correspondingly input of wave action from wind and its dissipation. This 
description is based on phenomenological parameterizations mainly [see [l3| . It gives very high dispersion of estimates 
of the terms [see fig. 1 in@] and a number of physically different compositions of S n i, Si n , Sdiss- Validity and physical 
correctness of the empirically based terms Si n and Sdiss are generally beyond critical consideration: quantitative 
aspects are dominating over obvious questions on their physical relevance. It is not the case of the theoretically based 
term S n i which explicit form provokes questions on limits of the asymptotic model of wave-wave interactions. 
For deep water waves four-wave resonant interactions are described by the collision integral 

S nl = Ttg 2 J\T 0123 \ 2 {N N 1 N 2 + N l N 2 N 3 -N N 1 N 2 - N Q NxN 3 ) 

(2) 

x S(k + ki - k 2 - k 3 )(5(u;o + u>i — u 2 — uj3)dk 1 dk 2 dk3 



Here we use convention [14j for normalization of spectral densities N(k) and kernels T0123. We also use subscripts 
i = ... 3 to denote dependence on argument [see [l5[ . Collection of cumbersome expressions for the interaction 
kernel To 123 can be found in Appendix of 0]. 



A. Homogeneity property of S„i and validity of asymptotic statistical description 

Homogeneity properties of S n i are of key importance for our consideration. First, wave frequencies and wavevectors 
in |j2} obey the power-law dispersion relation for linear deep water waves 

w(k) = y/tfk\. (3) 

Interaction kernel IT0123I 2 is also a homogeneous function of power 6, i.e. 

\T(Kk ,Kk 1 , K k 2 ,K,k 3 )\ 2 - K 6 |T(k ,k 1 ,k 2 ,k 3 )| 2 (4) 

Correspondingly, for the collision integral ([2]) that is cubic in wave action spectral density one has the well-known 
re-scaling property in terms of wavenumber dependence 

S nl [kk, vN k ] = k 19 / VS rii [k, A k ] (5) 

or in terms of dependence on frequency u> 

S nl [vu, vN^] = v n v s S n i[u, N u ] (6) 

One can see that the collision integral depends heavily on wave scales and amplitudes that may break validity of the 
asymptotic approach during long-time evolution of wave field. The basic assumption of the asymptotic approach is 
smallness of wave period T as compared with time scale of nonlinear interactions, i.e. 

T_ = _J_dN± = _J_ Sni <<; l (?) 
T n i WkA k dt WkA k 
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It can be checked that for 

iV k ~|k|- 9 / 2 or (8a, b) 



conditions (|5l6p are satisfied for any stretching parameter k and v. In other words, the asymptotic approach appears 
to be formally valid at any wave scale. Moreover, one can prove that this is true at any order of the asymptotic 
approach. For the collision integral asymptotic series 

oo 

Snl = J2 S % ] ( 9 ) 

n=4 

conditions similar to ([7]) remain valid for every term Sffi that represents resonant interaction of n waves [l|. Thus, the 
generalized Phillips spectrum ([SJ correspond to a special case when asymptotic description of wave-wave interaction 
can be formally extended over the whole range of wave scale. It gives a chance to propose a model where weak 
nonlinearity acts as a strong physical factor and is able to concur with other strong physical mechanisms, say, with 
dissipation due to wave breaking. 

It should be stressed that the found Phillips spectrum for deep water waves differs from the classic Kolmogorov- 
Zakharov solutions for direct and inverse cascading [H, 0] 

AT«(k) = C p P 1/3 g- 2/3 \k\- 4 (10) 

7V(2)( k ) = C g g 1 /3 g -l/2| k |-23/6 (n) 

(P, Q - energy and wave action flux, C p , C q — the corresponding Kolmogorov's constants). Collision integral S n i for 
solutions ([TUl ITT]) is plain zero and estimates of the kinetic equation validity requires special care. Following the 
simplest (but not trivial) way proposed recently 0, [Tij split nonlinear transfer term S n i into two parts - nonlinear 
forcing and nonlinear damping T^N^ (T^ - nonlinear damping rate) as follows 

Sni = Fk. — TkiVk 

Evidently, the relaxation rate Tk gives a physically correct estimate of time scale of nonlinear wave- wave interactions 
in the kinetic equation (p}. Thus, the KZ solution breaks when [see eq.17 in [t| 



that gives 



for direct cascade solution [5[ and 



T k LU ~ 47T ff |k| 9 iV2 ~ 1 

47rC p 2 5 - 1/3 ^ 2/3 |k b ,|^l 



4wC q Q 2 / 3 \k, 



he 



4/3 



for the inverse cascade @. The generalized Phillips' spectrum ([8]) being steeper than its KZ counterparts cannot be 
balanced 'on its own'. In contrast to KZ solutions (|10l lll|) it requires an additional term to provide a full balance in 
the right-hand side of the kinetic equation. 



B. Matching the dissipation function with nonlinear transfer term 



Matching a dissipation function can be done from different physical positions. First, one can rely completely on 
dimensional consideration and to construct a nonlinear dissipation function in the spirit of conventional models of 
wind- wave modeling as a local nonlinear function of spectral density AT(k) (N(ui)), frequency u> (or wavenumber |k| 
and gravity acceleration g as the only physical scale of deep water waves. To match homogeneity conditions ([7]) of 
the nonlinear transfer term for the Phillips spectrum (|8]) one has 

5 dlss [k,iV r ] = 7kN k = ^> k (g 1 ^ 2 N k )u J N k , (12) 



(13) 
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with non-dimensional constants Ck, C w and arbitrary functions $k or ^ of non-dimensional arguments. Such 
'conventional' approach implies that wave dissipation is determined by instantaneous state of wave field and, in a 
sense, represents the Phillips authentic idea of saturated spectra Q. 

An alternative way to introduce the nonlinear dissipation for the Phillips spectrum follows the Phillips idea 3] of an 
equilibrium of different physical mechanisms, in our case, wave dissipation and nonlinear transfer. It can be regarded 
as 'more physical': we assume that flux of energy (or wave action) due to nonlinear transfer is damped by nonlinear 
dissipation. Thus, one has a dynamical equilibrium provided by these two physical mechanisms. Such vision of wave 
balance can be formulated quite naturally in terms of wave turbulence where spectral flux of energy P plays a key 
role 

S diss [uJ, E u ] = Cph^iP^/g 2 ) - (14) 
where \l/ is arbitrary function of non-dimensional argument. 



III. THE PHILLIPS DISSIPATION FUNCTION AND PHILLIPS' SPECTRUM 

Consider the simplest model of balance of nonlinear transfer and wave dissipation (fT4|) that satisfy condition 
and provides the Phillips spectrum (jHJ formally in the whole range of wave scales 

dE(uj) dP(w) , , . 2 ,P(lo) , , 

dt = — ^~* (Pcj/5) 4r (15) 

The right-hand side of (1151) operates with formally non-local quantity of spectral flux of wave energy 



P(w) = -/ / Sni[u,6,E u ]dudB 

JO J -IT 



(16) 



In fact, this non-local value can be interpreted quite naturally using homogeneity properties ([5l|6]). For gravity waves 
P ~ E 3 (uj)u} 12 (see [10]) one has 

g 2 ~ g 6 — K 11 ) 

i.e. the non-dimensional argument of \& is proportional to wave steepness in power 6. Note, that we use differential 
frequency-dependent steepness ji w in (|17[) in contrast to integral steepness 

4 = ^ (18) 

9 

(uip spectral peak frequency) or cumulative steepness 

fil = g~ 2 / E{u)uj 4 du>. (19) 
Jo 

The steepness definition (|18p in terms of total energy E to t gives always finite values. Typically, for growing wind sea 
it is below 0.1 and is decaying with time or fetch. Two alternative definitions are useful for description high-frequency 
part of wave spectra. For weakly turbulent direct cascade solution (|10p the differential and cumulative steepnesses 
are infinitely growing functions [i w ~ /i c ~ ^fu. The Phillips spectrum E(uS) ~ u>~ 5 gives saturation of differential 
steepness /i w and logarithmic growth of cumulative one /i c . 

Our interpretation of function $(Pu 3 / g 2 ) equalizes two different models of wave dissipation: local in spectral 
density (|T2l [T3l and non-local flux-dependent one (fl4l . We fix this parallel below in notation 

Plo 3 

— 2~ = M- 
9 
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A. Stationary solutions of the generalized kinetic equation 



The stationary solution of (ITBl for an arbitrary dependence ^(/x) gives nothing but dependence of steepness fi = 
Pcu 3 /g 2 on wave frequency 

6dAi In ( HL) (20) 



o M (3 -*(//)) \U! 

A finite uj± in (|20|) is introduced from dimensional consideration only. The denominator of the integrand in (|20p 
must not be zero at finite frequency u) for the integrand convergence and, hence, for physically correct values of fi 
(strictly positive). Thus, function ^(n) has to tending to a constant value. Any cases but ^(/x) — > 3 give power-law 
dependence of /i on frequency w. As far as \i has a sense of wave steepness the case ^f(n) — > 3 is seen as the only 
physically relevant one. In this case the steepness for infinitely short waves is decaying with uj. 

Important results can be obtained for power-law dependence of function \E' on non-dimensional spectral flux Puj 3 /g 2 . 
Let 



<r 



Stationary solution of ([To]) obeys 



Thus, the dissipation function \& (Puj 3 / g 2 ) should be constant for the stationary solution for any exponent R in (|21[) . 
Note, that for R = in (|2ip the formal stationary solution gives an arbitrary exponent that is determined by coefficient 
C, i.e. 

P~LO~ C 



Mathematically, an arbitrary coefficient C can give an arbitrary exponent for frequency spectrum. But from physical 
viewpoint we should take care because at R = the key physical parameter - gravity acceleration g falls out of the 
model. To avoid 'the non-physical freedom' one should take a formal limit at R — > that fix exponent C = 3. 

For R ^ exponent (—3) and dissipation rate C are fixed by condition (|22|) . Surprisingly, the condition of 
stationarity does not depend on choice of power-like function ^, i.e. on dependence of dissipation on wave steepness. 

From homogeneity property (|17[) one gets immediately the Phillips spectrum 

E{u) = a Ph g 2 u- b (23) 
where the Phillips aph can be determined from the relationship between energy spectrum and spectral flux. 



B. Self-similar solutions for the generalized kinetic equation 



Stationary solutions cannot be considered as a final argument for the Phillips spectrum: solutions that evolve in 
time and space should be studied as routes to possible equilibria states of wave spectra. Fortunately, for R = the 
generalized kinetic equation (JT5J) has self-similar solutions and re-scaling property (| 1 T[) can be used in its full. 

Consider spatially homogeneous (the so-called duration-limited) solutions of (| 15[) in the form 

E - t p+ "U p (u)t q ) (24) 

Total energy of this solution grows as 

Etot - t p (25) 

Looking for similar homogeneity properties of terms in the right-hand side one gets from (|17[) the following conditions 
for exponents p, q 



9o- 1 



(26) 
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that is the same as for the family of self-similar solutions of the authentic conservative kinetic equation [see 0, for 
details] . The second condition single out special exponents for self-similar solutions of the generalized kinetic equation 

cm 

R(q -1) = (27) 

The exponent q = 1 is unrealistic as soon as it corresponds to very fast growth of energy (E to t ~ t ) while R = 
looks quite attractive for the model of wave dissipation. 
The advantage of the proposed dissipation function in (|15p 



Sdiss — —Cph — ^— (28) 

is two-fold. First, self-similar solutions are seen usually as physically relevant asymptotics of solutions with arbitrary 
initial and boundary conditions. The second feature of the dissipation function is more important for further study. 
While the nonlinear transfer and dissipation terms have similar properties of homogeneity they can be balanced 
without additional fitting parameters. Thus, choosing the dissipation function (f2"5|) we follow the Phillips idea 'like 
cures like'. 

Note, that Q realized this principle for all three terms in the right-hand side of the kinetic equation by fitting 
essentially nonlinear collision integral S n i and quasi-linear functions Si n , Sdiss to each other. The outcome of such 
composition was the well-known energy spectrum E(uS) ~ associated with a number of experimental findings [e.g. 
Il6l - ll8l | . Extremely high price for this sophisticated composition was a number of uncertainties in empirical parameters 
of functions of wave input and dissipation [see sect. 5 'Constraints on the empirical constants' inQ. 

Following the Phillips principle 'like cures like' for the case of self-similar evolution with dissipation function 
we avoid the Phillips problem: 

the spectral balance in our model does not require any fitting parameters. 



Universality of dissipation function for the Phillips spectrum E ~ cj 



-5 



Important analytical results can be acquired from analysis developed by [19| . Let us estimate collision integral S n i 
for power-like isotropic functions 



One can find an explicit expression 



E(u) = E [ — ) (29) 

UJQ 



5*? = SttV^o 11 (^) U ^ Fix) (30) 



where function F(x) is calculated by integration in parameter q - 'base vector' of a resonant quadruplet. Here we 
used superscript for collision integral in the kinetic equation (fTJ) written for energy. Function F(x) is shown in fig. Q] 
This function tends to infinity at x = 1 and x — 11/2 when the corresponding integral in parameter q diverges. The 
first case x = 1 gives divergence in high frequencies while x = 11/2 — in low frequencies. Two zeroes of function F[x) 
correspond to classic Kolmogorov-Zakharov solutions for direct (x = 4) and inverse (x = 11/3) cascades (see eqs. IT0| 

en). 

Similar analytical expression can be obtained for spectral flux of wave energy P 

'■^(s)"^ (31) 

For the inverse cascade case x = 11/3 the energy flux vanishes F(x) = - there is no but wave action flux to low 
frequencies. For x = 4 - the direct cascade F(x) tends to zero together with the denominator in (|3"Tj) that gives a 
finite constant of wave energy flux. 

Comparing (|30l and (|3"Tj) one can get a ratio that can be related with dissipation coefficient Cph in ([28]) for formally 
stationary solutions 



(32) 



FIG. 1: Function F(x). In the right panel the plot is zoomed in order to show its behavior near points corresponding to the 
Kolmogorov-Zakharov solutions x — 11/3 and x = 4. 



For the Phillips spectrum with x — 5 this ratio gives Cph = 3, i.e. exactly the coefficient for an arbitrary dependence 
of dissipation function on non-dimensional argument Pw 3 /g 2 (see!22p. 

Below we present results of simulations where dissipation will be expressed by 'conventional' dependencies on 
spectral densities (|12ll3p or differential wave steepness ([TT]). For 

Sdiss[u,E(u))] = C 'pumps n%u)E{ui) (33) 

this dissipation function balances the explicit collision integral (|30[) when 

Cp h unps = 8tt 2 F(5) « 2.03 (34) 

Here we used an approximate expression for F(x) at x — ► 11/2 

, . 0.0129 
F(x) = TT7 ^, x-11/2 

Our estimate assumes the balance of nonlinear transfer and dissipation in the whole frequency range that is not the 
case of wind waves where strong inherently nonlinear dissipation occurs for relatively short waves as a result of wave 
breaking. More realistic dissipation function with a cutoff fixed frequency uidiss 

Sdiss = CdissH%uE{bj)@{bJ - uJdiss) (35) 

takes this phenomenon into account in the simplest way. Here Q(x — x ) is Heaviside function. This form of 
the dissipation function will be used below in our simulations. The theoretical estimate of the Phillips dissipation 
coefficient Cphunps (1341) gives a good reference in our numerical study. 

IV. NUMERICAL RESULTS 

Numerical simulation of the kinetic equation ([1} has been carried out in order to check the simple theory presented 
above. The WRT algorithm [2(J[2l[ was used for calculation collision integral. The dissipation function was calculated 
in accordance with definition (|33|) . i.e. in terms of 'observable' parameters of wave steepness. In this paper we present 
the very hrst results of simulations in the simplest setup. We studied duration limited growth (spatially homogeneous) 
of initially isotropic spectrum. The initial spectrum is given by JONSWAP-like formula that correspond to frequency 
spectrum of energy 

^ = ^ 4eXP ("l(S)") 7G (36) 

where 
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FIG. 2: The kinetic equation solutions for reference case of swell: top-left - frequency spectra of energy, hard straight line - 
E(w) ~ u)~ 4 , dashed - to~ 5 ; top-right - energy flux (curves), straight line corresponds to law P ~ uj~ 3 ; bottom-left - cumulative 
(lines) and integral wave steepness (symbols at u)/lj p = 1); bottom-right - differential steepness defined in terms of spectral 
flux (bold line) and parameter fj, w . Time in seconds is given in legend. 



with default peakedness parameters 7 = 3.3, a p . The initial wavelength corresponds to rather long swell A « 240 
m. Minimal frequency in simulations f m in = 0.02Hz, maximal - f m ax = 2Hz, peak frequency f p = 0.079Hz, the 
dissipation cutoff is set at four times higher frequency fdiss = 0.319Hz. With logarithmic grid each frequency domain 
(low bound / m i„-peak, peak-dissipation cutoff) contains approximately the same number of grid points. In absence 
of wave input slow spectral downshift does not change essentially this division. 

The reference case of 'pure swell' is presented in fig. [5] As it was demonstrated in Q the solution evolves to a 
self-similar shape. In our case we see a transition to a self-similar behavior rather than the self-similarity itself because 
of relatively short time of evolution and realistic low wave steepness. Energy spectra (top fig. [2]) show downshifting 
and a peakedness growth. The energy spectral flux tends to constant magnitudes in a wide range of non-dimensional 
frequencies (3 < lo/uj p < 20). Integral wave steepness [i p (circles in bottom fig. H]) is going down with time as well 
as curves of cumulative steepness /i c which is growing function of frequency. The differential steepness shows similar 
behavior - steep growth with frequency. 

Fig. [3] present results for dissipation function with coefficient Cdiss = 1-22 = Q-QCphaups and dissipation cutoff 
frequency todiss = 2rad _1 . This coefficient is almost two times less than our estimate of the Phillips coefficient 
Cphiiii ps — 2.03 defined for the Phillips exponent (—5). One can see appearance of essentially steeper spectrum tail 
and dramatic reduction of all the steepness parameters defined above while for cumulative steepness fi c one can see 
a sort of saturation. At the same time, the differential steepness fi w shows a gradual growth with frequency. In 
bottom-right fig. [3] we give curves in accordance with two possible definition of dissipation, first, in terms of spectral 
flux (I14[) . second, as a function of wave steepness fj, w (1331) and conventional presentation of dissipation in wave models 
([T2l IT3l) . Higher dissipation rate Cph — 1.9 - ~ 5% lower than the Phillips dissipation coefficient Cphiiu ps = 2.03 
shows very strong tendency to Phillips spectrum E(iS) ~ uj~ 5 . The bottom-right panel fixes two key points: first, 
the differential steepness is tending to saturation, second, two alternative definitions of the parameter, in terms of 
spectral flux (bold lines) and as function of wave amplitude (thin line), appear to be close to each other. 

One more example illustrate the effect of dissipation cutoff. In fig. [5] all the parameters of simulation are the same 
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FIG. 3: The kinetic equation solutions for dissipation function (|33l) with Cph = 1-22 and dissipation cutoff tvdiss = 2rad _1 : top- 
left - frequency spectra of energy, hard straight line - E(ui) ~ u>~ 4 , dashed - uj~ 5 ; top-right - energy flux (curves), straight line 
corresponds to law P ~ bottom-left - cumulative (lines) and integral wave steepness (symbols at uj/lu p — 1); bottom-right 
- differential steepness defined in terms of spectral flux (bold line) and parameter fi w . Time in seconds is given in legend. 



as in previous case but dissipation cutoff is set at Udiss = 3rad . One can see the same effect: very strong tendency 
of solution to Phillips' asymptotics and saturation of all the steepness parameters (fx p , (x c , fi w ) 



V. CONCLUSIONS AND DISCUSSION 



In this work we proposed to reanimate the Phillips idea of wave spectra as a result of a balance of different physical 
mechanisms. We combine effects of nonlinear transfer due to four- wave resonant interactions of weakly nonlinear deep 
water waves and dissipation dealing with wave breaking. These two effects can be co-existing in quite good symbiosis 
for the famous Phillips spectrum E(uj) ~ (JSJ) : just this spectral shape makes asymptotic kinetic equation to be 
formally valid in the whole range of wave scales. Thus, we have a good example of the generalized Phillips spectrum 
in this case. 

There is one more point of interest in competition of dissipation and nonlinearity for the Phillips spectrum. Sta- 
tionarity of this spectrum requires certain magnitude of dissipation coefficient - the Phillips coefficient. This is in 
contrast to widely accepted models of wave dissipation where dissipation function can be arbitrary: the spectral shap- 
ing does not depend on the dissipation magnitude. We gave analytical estimate of the Phillips dissipation coefficient 
C phiiups ~ 2.03. 

Our simulation within the simplest setup of slowly evolving isotropic wave field justified our theoretical results. 
Shape of wave spectra are shown to be sensitive to magnitude of wave dissipation. For values close to the found 
theoretical value Cphuu P s ~ 2.03 the spectra are tending to the Phillips classic tail w~ 5 and appeared to be long- 
lived. Deviation from this value give other spectral slopes. Variety of spectral slope is well-known from experimental 
studies. Collection of spectral slopes derived from buoy data have been presented, say, by Paul Liu long time ago 
[ill [l2l ]. His fig. 8 in [l2| shows a wide range of the exponents: 90% of them fall into range 3 — 7 with a maximum 
between 4 and 5. This high dispersion has been attributed to 'an equilibrium range' of wave spectra with no clear 
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FIG. 4: The kinetic equation solutions for dissipation function (|33|l with Cph = 1.9 and dissipation cutoff Udiss = 2rad~ : top- 
left - frequency spectra of energy, hard straight line - E(ui) ~ u>~ 4 , dashed - lu~ 5 ; top-right - energy flux (curves), straight line 
corresponds to law P ~ bottom-left - cumulative (lines) and integral wave steepness (symbols at uj/lu p — 1); bottom-right 
- differential steepness defined in terms of spectral flux (bold line) and parameter fi w . Time in seconds is given in legend. 



idea what mechanisms are responsible for this equilibrium and why the Phillips spectrum cj~ 5 is emphasized in this 
experimental collection. 

Our very first numerical experiments with the new dissipation function showed robustness of the effect of nonlinear 
dissipation. Its effect leads to saturation of wave steepness both in integral (fj, p ) or differential (/i w ) quantities. This 
feature is of great importance for further simulation of wave spectra with both effects of wave input and dissipation 
taken into account. High nonlinearity of the new dissipation function suppresses effect of wave pumping in high 
frequency completely and, thus, the classic Phillips spectrum continues to play important role in general case. 

The research was conducted under grant of the Government of the Russian Federation 11.G34.31.0035 and ONR 
grant N000141010991. This support is gratefully acknowledged. 
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FIG. 5: The kinetic equation solutions for dissipation function (|33|l with Cph = 1.9 and dissipation cutoff Udiss = 3rad~ : top- 
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